Suppression of superconductivity by Neel-type magnetic fluctuations in the iron 

pnictides 
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Motivated by recent experimental detection of Neel-type ((7r,7r)) magnetic fluctuations in some 
iron pnictides, we study the impact of competing (n, n) and (n, 0) spin fluctuations on the super- 
conductivity of these materials. We show that even short-range, weak Neel fluctuations strongly 
suppress the s + ~ state, with the main effect arising from a repulsive contribution to the s" 1 pairing 
interaction, complemented by low frequency inelastic scattering. Further increasing the strength of 
the Neel fluctuations leads to a low-T c d-wave state, with a possible intermediate s + id phase. The 
results suggest that the absence of superconductivity in a series of hole-doped pnictides is due to the 
combination of short-range Neel fluctuations and pair-breaking impurity scattering, and also that 
T c of the pnictides could be further increased if residual (n, n) fluctuations could be eliminated. 



o 
o 

u . 
=5 

CZ3 ■ 

a- 



-a 

o 
o 



> 

(N 

cn 

O 

(N 



X 



The proximity of the superconducting state (SC) to a 
"stripe" spin-density wave instability (SDW) in the phase 
diagrams of the recently discovered iron-based supercon- 
ductors [l| (FeSC) prompted the proposal that SDW spin 
fluctuations provide the pairing mechanism 0|. Indeed, 
the Fermi surface of many iron pnictides consists of sep- 
arate electron and hole pockets with the electron pockets 
displaced from the central hole pockets by the SDW or- 
dering vector Qsdw = (tt>0) / (0, 7r) (see Fig. [Tj). In this 
situation, even weak SDW fluctuations may overcome a 



in the magnetic susceptibility at QnccI, which is however 
weaker than the peak at Qsdw [Z|- 

Recent experiments also evidence the importance of 
Neel order and fluctuations in some FeSCs. Neu- 
tron scattering measurements @ revealed that even at 
small x, Ba (Fei_ x Mn x ) 2 As2 exhibits spin fluctuations 
peaked at QnccI, m addition to the SDW fluctuations 
peaked at Qsdw- NMR measurements Q confirmed that 
these Neel fluctuations couple to the conduction elec- 
trons. Because the entire family of "in-plane" hole-doped 



strong on-site repulsion giving rise to an s+ SC state, Ba (Fei_ x M x ) 2 As 2 compounds (M = Mn, Cr, Mo) 
in which the gap function has one sign on the electron 
pockets and another sign on the hole pockets Q. 

As can be seen from Fig. [IJ the two electron pockets 
are connected by the momentum QnccI = (f, tt) suggest- 
ing that Neel-type magnetic fluctuations may also be im- 
portant These fluctuations favor a d-wave SC state 
in which the gap function has opposite sign in the two 
electron pockets. On the theory side, first-principle and 
Hartree-Fock calculations find that the Neel state is lo- 
cally stable, but with a higher energy than the SDW state 
[H, 0|) while random phase approximation (RPA) calcu- 
lations performed in the paramagnetic phase find a peak 
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Figure 1: (left panel) Schematic Fermi surface configuration 
in the 1-Fe Brillouin zone, with two central hole pockets and 
two electron pockets, (right panel) Self-energy diagrams of 
the Eliashberg equations: normal component (upper panel) 
and anomalous component (lower panel). 



displays SDW order at x = and Neel order at x = 1, 
we expect that competing Neel and SDW fluctuations 
will be found across the whole material family, making 
this series a promising system to tune these competing 
spin fluctuations. While superconductin g sta tes have not 
been reported in these materials to date [Tl| (in contrast 
to the electron-doped counterparts M — Co, Ni, Rh, Pt, 
Cu, where SC is always observed (l2|) higher pressure ex- 
periments on samples with lower residual resistivity may 
reveal SC. 

There is also indirect evidence for Neel fluctuations 
in the extremely electron-doped A y Fc2- x Sc2 compounds 
[l3j . In these materials, the near absence of Fermi sur- 
face pockets in the center of the Brillouin zone suggests 
that SDW fluctuations and the s" 1 state are disfavored, 
while the square-like shape of the electron pockets is ex- 
pected to enhance the (7r,7r) fluctuations [1J|- Chemical 
substitution on the A site or application of pressure (l5| , 
however, can create a small pocket in the center of the 
Brillouin zone, which could support (tt, 0) fluctuations 
and s" 1 superconductivity. 

The effect of competing spin fluctuations on FeSCs is 
thus of experimental and theoretical interest. In this pa- 
per, we address the problem via a multi-band Eliashberg 
approach fltl [l7| in which the spin fluctuations are mod- 
eled as bosons with spectral density determined experi- 
mentally or chosen phenomenologically, while their effect 
on the electrons is determined from the one-loop approx- 
imation to the self energy (see diagrams in Fig. [T]). This 
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approximation has been extensively employed in studies 
of cuprates [l8l - l20| , ferromagnetic SC [2l| , and pnictides 
(2^ |. Although some complications arise for two dimen- 
sional electrons coupled to two dimensional spin fluctu- 
ations, it is generally accepted that the one-loop calcu- 
lation gives a reasonable representation of the physics 
when the typical boson frequency is small compared to 
the Fermi energy (20| . 

Our calculation goes beyond previous work [22J by in- 
corporating both SDW and Neel fluctuations and includ- 
ing the Coulomb pseudo-potential, as well as by using 
the experimentally determined spin fluctuation spectrum 
instead of the single-pole approximation employed previ- 
ously. We find that the Coulomb pseudo-potential has 
only a weak effect on the dominant state but that 
even weak, short-range Neel fluctuations strongly sup- 
press the transition temperature y^~ wavc . if sufficiently 
strong, the Neel fluctuations may induce a d-wave state, 
but the transition temperature is found to be much lower 
than the optimal T c for the s H state. The transition be- 
tween s" 1 and d-SC may either occur via an intermediate 
time reversal symmetry-breaking s + id state (23|, [24j or, 
if the impurity scattering is stronger, via an intermediate 
non-SC state separating the two regions (see Fig. [2]). We 
further show that the low temperature state existing 
in the presence of combined SDW and Neel fluctuations 
is easily suppressed by impurities or, presumably, by a 
competing ordered state, such as the SDW [2^, [26| . 

To gain insight into the results, we use functional 
derivative methods pioneered by Bergmann and Rainer 
[27l |28| . We find that the strong suppression of the s H 
state comes mostly from a repulsive s" 1 pairing inter- 
action induced by the Neel fluctuations, although pair- 
breaking inelastic scattering plays some role. Finally, we 
discuss the implications of our results not only to the SC 
of the in-plane hole-doped pnictides, but also to the value 
of T c in the FeSCs in general. 

Our model consists of a two-dimensional Fermi surface 
with two hole pockets (r, density of states N-p) at the 
center of the Brillouin zone and two electron pockets (X 
and Y, density of states Nx) displaced from the center by 
the momenta (71", 0) and (0, n) (Fig. [T]) [2{|. For simplic- 
ity, hereafter we assume that these two hole pockets are 
degenerate - our results do not depend on this simplifica- 
tion. Following Ref. [H, we set r = N x /N r = 0.65. The 
electrons are coupled to two types of low-energy bosonic 
excitations, namely, SDW spin fluctuations peaked at 
(it, 0) / (0, 7r) and Ne el sp in fluctuations peaked at (ir, it). 
Experiment (Refs.pl. l3l|) indicates that in the paramag- 
netic phase these excitations are described by diffusive 
dynamic susceptibilities: 
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Here, q is the momentum deviation from the ordering 
vector Qi (all lengths are in units of the lattice parame- 
ter a) and f2„ is the bosonic Matsubara frequency. The 
quantity that actually enters the Eliashberg equations 
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Figure 2: Transition temperatures T c of the s-wave (red/light 
curve) and d-wave (blue/heavy curve) states as function of 
the Neel magnetic correlation length $Neel : for r = 0.65, 
7ncci/7sdw = 0.33, ANeei/AsDw = 2, and fi* — 0.8. T c ,o « 
O.I7SDW is the transition temperature for £nc C i = 0. The 
shaded area denotes the regime where the two states have 
similar transition temperatures and a possible s + id state 
may occur. The dashed lines show the behavior of the system 
in the presence of impurity scattering. The inset shows the 
frequency dependence of the spectral function ^Neei (w) of the 
Neel fluctuations for different values of ^Neei- 



is the spectral function integrated over the momentum 
component g|| parallel to the Fermi surface and evalu- 
ated at q ± = 0, i.e. A Ncc \(cu) = f dq\\lm.XNeel (q\\,w) ■ 
This spectral function gives rise to the Matsubara-axis 

interaction a^ 1 ' (f2 n ) 



£i/\Jl + |fi„| 7i 1 Q which enters 
the Eliashberg equations as described below. The spin 
fluctuations in each momentum channel i are described 
by two parameters: the Landau damping 7$, which sets 
the energy scale of each fluctuation spectrum, and the 
correlation length £j, which sets both the strength and 
the spatial/temporal correlations of the spin fluctuations. 
We will tune the spectrum by varying 

Because the Landau damping originates from the low- 
energy decay of the spin excitations into electron-hole 
pairs, 7, is determined by the electron-boson coupling 
constant gi and the densities of states. Within a stan- 
dard spin fluctuation model there is one adjustable pa- 
rameter, which we choose to be gsDW an< i which we set 
to yield j^~ wave ~ 30 K. The remaining quantities are 
determined from the experimentally measured 7^ and the 
ratio of densities of states r (see the Supplementary Ma- 
terial for details). Following the experimental results of 
Ref. @, we use 7n co i/7sdw ~ 0.33 with 7sdw ~ 25 
meV; the value of (?Neei follows from the relationship 
between 7n ci/7sdw and .9n oc i/.9sdw- Finally, we set 
£sdw = 5a throughout our calculations, varying the cor- 
relation length of the Neel fluctuations £ncc1- Our results 
do not change significantly for smaller values of £sdw- 
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To obtain the transition temperatures in the s and 
d-wave channels we linearize the Eliashberg equations 
in the superconducting quantities and solve the result- 
ing equations for the anomalous component W an and 
the normal component Z ajl = ImS^ n /iuj n of the self- 
energy (the real part of T, N just renormalizes the band 
dispersions, possibly differently for different pockets [Trl 
[3H3)- These quantities are averaged over each Fermi 
pocket becoming functions only of the Fermi pocket la- 
bel a and the fermionic Matsubara frequency u) n = 
(2n + 1) 7rT. With the aid of the auxiliary "gap func- 

W r .„ i a _ W { x/Y), n 



tions" A 



and A 



Zr, n \u n \VNx W*J>« — Zx^u^VNr' 

the linearized gap equation is expressed as a matrix equa- 
tion in Matsubara (indices n, m) and band (indices a, /?) 
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Here we have introduced the matrix elements com- 



ing from the bosonic modes a 
1 corresponds to SDW and i 



(i) 
nm 



= a l (uj„ - co m ) (i = 
2, to Neel fluctua- 



tions) and the dimensionless coupling constants Asdw = 
2 .9sDW\/^riVx, AncoI = S^eei-^- T is tne temperature. 
We also introduce an upper frequency cutoff A = 87SDW; 
corresponding to the energy scale of the bottom/top of 
the electron/hole bands, and we assume that fj,* is a bare 
Coulomb interaction renormalized in the standard way 
by higher energy processes, r -1 is the scattering rate 
associated with non-magnetic point impurities and the 
Z a n functions are obtained analytically (see Supplemen- 
tary Material). The Coulomb pseudo-potential favors so- 
lutions with J2 a N a A a = 0. 

Reflecting the tetragonal symmetry of the system, this 
matrix equation supports two different types of solu- 
tion: the s-wave state Ax,n = Ay j7l , with either ,s ++ 
(Ar.n oc Ax,n) or s_l (Ar.n oc —Ax, n ) structure, and 
the d-wave state Ax. n — -A^„. The solution in a given 
symmetry channel is obtained when the largest eigen- 
value 77 of the matrix @ vanishes. Since our calculations 
never yield an s ++ state, hereafter we use the terms s- 
wave and to refer to the same state. Due to limi- 
tations of the size of the matrices that can be diagonal- 
ized, and since the matrix size scales as A/T, we resolve 
T c > 10 _3 7sdw- Hereafter, we set Asdw = 0.4 and the 
Coulomb pseudo-potential fj,* = 0.8, which gives, in the 
absence of competing Neel fluctuations, j^~ wavc ~ 30K 
and implies AnccI = 0.8. 

Fig. [2] shows our principal results: the dependence of 
the SC transition temperature T c on the strength of Neel 
fluctuations (parametrized by the Neel correlation length 
^Neei)- The light solid line (red online) shows the tran- 
sition temperature ys-wavc f Qr ^ ne s h channel in the 
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Figure 3: Averaged s + ~ gap function ^2 N a A a across the 
different pockets at j 1 =- wavc as function of Matsubara fre- 
quency u) n (in units of 7sdw), for fx* = (green/dashed curve) 
and // = 0.8 (red/solid curve). The inset shows T c S Q Wavc (in 
units of 7sdw) as function of fi* . Although here we used 
?Neei = 0, a similar behavior holds for ^Neei 7^ 0. 



absence of impurity scattering. One sees that even weak, 
short-range fluctuations strongly suppress SC, but 
that once T^~ waMC has been substantially reduced, the 
additional suppression effect caused by further increas- 
ing £ncc1 is small. Sufficiently strong Neel correlations 
produce a d-wave solution (heavier solid line, blue on- 
line) with j'd- wavo that eventually becomes larger than 
ys-wavc arwa y S remains small compared to the max- 
imum ys-wavc^ j n Qur Xiriearized theory the transition 
between s-wave and d-wave superconductors appears as 
a discontinuous ch ang e in the nature of the state, but the 
considerations of [23j suggest that the nonlinear terms 
not included here change the physics so that the s and 
d states are separated by a range of s + id supercon- 
ductivity (shaded area). The dashed lines show the be- 
havior in the presence of impurity scattering, which is 
pair-breaking for both and d-wave superconductiv- 
ity. Sufficiently strong impurity scattering can discon- 
nect the two SC states, leaving an intermediate non-SC 
regime. 

We also analyze the impact of the Coulomb pseudo- 
potential fi* on the s H state - the d-wave state avoids 
the Coulomb repulsion. Fig. |3]shows the pocket-averaged 
gap function ^ Q iV Q A Q both in the presence and in 
the absence of fi*. To avoid the local repulsion, the av- 
eraged order parameter changes sign at a non-zero Mat- 
subara frequency, although the sign of each individual 
gap does not necessarily change. This is the multi-band 
analogue of the response of a single-band s-wave super- 
conductor to the local repulsion. For all values of £ncc1 
we have studied, neither A n= o nor ys-wavc (shown in the 
inset) are substantially altered by //*, in agreement with 
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Figure 4: Functional derivative T^uST^'™*** /6A Ne ei (w) as 
function of frequency w (in units of 7sdw) for the cases 
^Ncoi = (red/solid line), £n C c1 = 0.1a (green/dashed line), 
and £nco1 = 0.5a (blue/dotted dashed line). 



the weak-coupling analysis of Ref. [33] . 

We now turn to the physics of the decrease of T*~ wav ° 
caused by Neel fluctuations. Increasing ^Neei increases 
the spin fluctuation intensity and changes its functional 
form(see the inset of Fig. To understand in more 
detail the effects on pairing, we use the func- 

tional derivative techniques introduced by Bergmann and 
Rainer [2?], HH to analyze how different frequency regions 
of the Neel spectral function affect T^ wavc . In Fig. H 
we plot the functional derivative 
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for different values of ^Neei- Previous work [271 . l28fl has 
shown that the low frequency regime captures the pair- 
breaking effects of inelastic scattering while the high fre- 
quency regime expresses the changes to the pairing in- 
teraction. The larger magnitude of <51ogT c /(5AN CC i at 
high frequencies shows that in the pnictides the domi- 
nant effect of the Neel fluctuations is to provide a neg- 
ative contribution to the pairing interaction, with 
the extra pair-breaking effect of the induced low fre- 
quency inelastic scattering being non-zero but less impor- 
tant. Because Fig. @] shows that the logarithmic deriva- 
tive of T c is a slow function of ^Neei we conclude that 
the initial steep drop and subsequent flattening of the T c 
curve shown in Fig. [2] is due in large part to the varia- 
tion of T c itself. Additionally, as ^Neei is increased the 
Neel fluctuation spectrum shifts to lower frequencies (see 
the inset of Fig. [2]), where the pair-breaking is less effec- 
tive. However, additional physics is also at play. Con- 
sider, for instance, the weak-coupling limit of two effec- 
tive competing pairing interactions A s and A^. We obtain 



V(v / Af+Al-A d ) 
dlogT c 



and 



A? + ^-A dv /AjTAl) 



(4) 



which is larger in magnitude for larger A^, implying an 
opposite ordering of the curves to that seen in Fig. [4] 
The difference between our results and this simple weak- 
coupling estimate arises because the gap function self- 
consistently adjusts to the pairing potential, so that for 
larger £ncc1 the gap function decreases more rapidly with 
frequency, thereby minimizing the depairing effects of the 
Neel fluctuations (see Supplementary Material). 

Our results offer a possible explanation for the puz- 
zling behavior of the hole-doped Ba (Fei_ x M x ) 2 As2 se- 
ries (M = Mn, Cr, Mo) [loj], which, in contrast to its 
electron-doped counterpart (M = Co, Ni, Rh, Pt, Cu) 
fl^ . does not display SC. The short-range Neel fluctu- 
ations induced by the dopants, which were observed ex- 
perimentally for low concentrations of M = Mn [f| Q , 
suppress the s" 1 state without giving rise to a high- 
temperature d-wave state. This low-T c state, in 
turn, can be easily suppressed, for example by impurity 
scattering or by another competing ordered state, such 
as the SDW state j2f| observed at low x. We suggest 
that improving the purity of the samples and applying 
pressure to suppress the SDW state may reveal either a 
weakened state or perhaps a low T c d-wave state. 
Similarly, in the extremely electron-doped A y Fe2_ x Se2 
systems [13l . l3~il , Hf| , where for small y the hole pocket 
is generally absent and d-wave superconductivity is dis- 
cussed, adding holes by changing A [l4j] or by applying 
pressure (l5| should produce the reverse competition. 

More generally, we expect that most of the FeSC fam- 
ilies will display Neel-type fluctuations due to the ex- 
istence of two matching electron pockets separated by 
QNeei- Our findings show that, even when they are weak, 
these (7r, 7t) fluctuations will strongly suppress J^ _wavc . 
This suggests that one possible route to optimizing T c is 
to minimize these fluctuations, for example by making 
the two electron pockets unequal via a tetragonal sym- 
metry breaking [36]. Interestingly, torque magnetometry 
measurements found such a tetragonal symmetr y b reak- 
ing above T c in some optimally doped FeSCs j37(. In 
Ref. [HI , it was also observed that a small strain applied 
along the orthorhombic axis can enhance T c . 

In summary, our results open a new route to explore 
unconventional superconductivity in multi-band systems 
by controlling competing spin fluctuations. In particu- 
lar, Neel fluctuations have a strong effect on the s H 
state of the FeSCs, rapidly reducing T c and potentially 
driving a transition from s-wave to d-wave SC (Fig. [2]). 
Depending on the strength of impurity scattering, more 
exotic states can emerge, such as the s + id state [23| . 
Note that an s + id state was also suggested to emerge 
due to other mechanisms, such as extreme doping and 
Fe-As hybridization [24]]. Notice also that, in our model, 
since there are two SC instabilities in different symme- 
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try channels, the lower T c solution, while not present in 
the ground state, will give rise to a collective excitation 
which can in principle be detected by Raman scattering 
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Supplementary material for "Suppression of superconductivity by Neel- 
type magnetic fluctuations in the iron pnictides" 



I. FORMULATION AND SOLUTION OF ELIASHBERG EQUATIONS 



A. Formulation of Equations 



The low-energy action describing the coupling between the electrons and the SDW and Neel fluctuations is conve- 
niently expressed in terms of the Nambu operator = ( cj, 

where cr.k<r, cr',k(D cx^a, and cy^o- correspond to operators on the T and T' hole pockets, on the X electron pocket 
at (it, 0) , and on the Y electron pocket at (0,7r), respectively. We have: 



p Op Op 

s= hf{ (e k - ujJ) ^ + E / ( fc ) • + E si / *i,-k-« • (*lf>i%) (si) 
Jk i=1 j k i=1 Jk, q \ ' 

where k = (k, u>„) refers to both momentum and fermionic Matsubara frequency, $i t k denotes the collective bosonic 
fields associated with the SDW (i = 1,...,4) and Neel (i = 5) fluctuations, and %% (&) refers to the corresponding 
dynamic magnetic susceptibilities. Our indices are defined such that X(r,r') (t, 0) corresponds to i = 1, 2, X(r.r') (0; tt) 
to i = 3,4, and x C 71 ") 7r ) to i = 5. The coupling constants satisfy gi=\..A = .9sdw an d 55 = (?Neei- We also have 
the band dispersions = diag 4 (er,k, £r',kj £x,k ; £y,k) ® T3, where Tj are Pauli matrices in Nambu space. For a 
spin-rotationally invariant system, and for the case of singlet pairing, we can focus on the z-axis projection of p i 
given by: 
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The Eliashberg equations are obtained by calculating the one-loop self-energy 

Jq 

with G7/ 1 = Gq k - E fc and Gq fc = - £ k . 



(S3) 



B. Reformulation of Equations 

To solve the self-consistent system of equations, Eq. lS3[ we rearrange them into a form more convenient for numerical 
solution. We write the self energy as 



S fe = iu n I 1 - Z k ) <g> r + Wk ® n + Cfe ® 7-3 



(S4) 



where Z fc = diag 4 (Z r , fe , Zr'.k, #x,k, 2y;fc) and Cfc = diag 4 (Cr,fc, Cr'.fe, Cx.fc, Cr.fc) are the imaginary and real parts of 
the normal component, respectively, and Wk = diag 4 (Wr.k, Wr',k> Wx,ki Wyk) is the anomalous component of the 
self-energy. As we discussed in the main text, the real part Ca.k renormalizes each a band dispersion (possibly in 
different ways Q), and will not be discussed here. 

Hereafter we will consider two degenerate hole pockets (Nr = Njv and xr (Qt) = XT' (Q<))j implying Wr = Wr 1 
and Zr = Z^i. Using the tetragonal symmetry of the system, we have x (0, 7r) = xC^jO) and = Zy, reducing 
the number of self-consistent equations to five. We integrate over the momentum component q± perpendicular to 
the Fermi surface, using the fact that the electronic propagator is more sharply peaked at the Fermi level than the 
bosonic propagator Q. Next, we linearize the equations by keeping the leading terms in order W a ,n and average 
the gaps along each Fermi pocket. We also include the impurity scattering and the Coulomb pseudo-potential in the 
standard way, obtaining, for the normal part: 



%9sdw n xT ^2 s S n ( 2m + x ) / d 1\\ Xsdw (q\\ , u n - u m ) 



(Zr,n — 1) w n = 



(Zjc,n-l)w„ = T^sgn(27«+1) 

m 

+9 2 u ^n YjN a + uf mp sgn (uj n ) N a 

a a 

and for the anomalous part: 



2 3sdw^t / dq\\Xs-Dvr(q\\,u n -u m ) +g$ leel Nx / dq\\XNeei (q\\, u n - u, 



(S5) 



11; 2 AT rpST (W X , m + W Y , m ) f ( v 

Vvr,n = -g S Dw N xl > ^ ; ; — / dq\\XsT>w{q\\,Un-Um) 

~~ £>X,m \u m \ J 



W a , m 2 ^ W a , n 



Za,m | w m| ^Q,n |^n| 

W(X/Y-),n = -2.9sDW Ar rT V „ Wr ; m 1 / dqu xsdw (gii , w„ - u m ) - g^ ccl N x T V ),m / dtfi XNeel , w n - w TO ) 

— ^r.m Wm J — AX.m Wm\ J 

m ' 1 ' m ' 1 1 



Here, g u is the coupling to the Coulomb repulsion and uf mp is the averaged local impurity potential. The Coulomb 
repulsion renormalizes all bare interactions, which become gf —> g\j (l + g^N), where N = 2Nr + 2Nx is the total 
density of states. 

Using the diffusive expression for the spin susceptibility, Eq. (1) of the main text, yields: 
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= 7rf. 
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i + |n»l7i e 
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For convenience, we absorb the coefficient ir into the coupiing constants gf and introduce the density of states 



ratio r = Nx/Nr, defining the Coulomb pseudo-potential /i* = g^N/ 



and the coupling constants A 



SDW 



2gl BVf y/N T N x and A 



Ncol 



5Ncc1 



1 + g^N), the scattering rate r 1 = 2Nxuf mp , 
Nx ■ An important note about the form of the 



magnetic susceptibility: while its static part comes from high-energy modes not considered in our model, the Landau 
damping is a direct result of the coupling between the paramagnons and the fermions, as given by Eq. (|S1|) . Thus, the 
parameter 7$ contains information about the coupling constant g^. By evaluating the bosonic self-energy to one- loop, 
we obtain 7sd W = K sDw5sDW-^f -^r and T^eei = K Neci5Ncci-^3o where Ki are dimensionless parameters presumably of 
similar orders of magnitude. This puts constraints on the ratio between the effective couplings Asdw — 2g| DW \/Ar^Vx 

and AnccI = SnccI-^^ > which is expressed as 



7SDW 
2«SDW 7Neel 



C. Solution of Equations 



Substituting the definitions given above into Eq. (|S5|) and evaluating the sums over Matsubara frequencies yields 



T 



T 



= (2n + 1) + \A*^SDW?SDW'5'sDW,ri + -~-sgn (uj n ) ( 1 + 



with the auxiliary functions: 



, A SDW ^ DW 5sDw,n . , , a \ . r 1 t \ ( i , 1 

(2n + 1) + ( -= — + ANeelSNeel^Neel.n I + "17r s g n ( w «) I 1 + ~ 



(S8) 



2 sgn {n) 



7« 



Here, Hw 
After 



sgn (n) + 1 7i 



2sgn(ra) + 1 , n = 0, -1 



+ sgn (n) , 0, -1 



2 ' 

defining 



is the Hurwitz zeta function for which efficient numerical evaluations exist, 
the auxiliary "gap functions" Ar, n = Wr. n / (V^x^r.n l w n|) 



(S9) 



and A 



(X/Y),n 



W(x/Y).n/ (VNrZx.n l w n|) and using the solutions for Z given above, it is straightforward to write down the 
gap equations (|S6[) as a matrix equation in Matsubara and band spaces, yielding Eq. (2) of the main text. For 
numerical computations, it is convenient to use the tetragonal symmetry of the system and split the matrix equation 
in two: one for the s-wave case Ax.n = Ay jn and another one for the d-wave case Ax, n = — Ay jT1 , yielding two 
different kernels K s (2x2 matrix) and Kd (lxl matrix). We use standard routines to obtain the leading eigenvalue 
and corresponding eigenvector of the matrix; the transition temperature is the temperature at which the leading 
eigenvalue crosses zero. 



II. FORM OF THE REAL AXIS SPECTRAL FUNCTION 



The spectral function for Neel fluctuations that enters the Eliashberg equations is 

^Nccl M = / 'dgylrnxNeei (?|h w ) 

Evaluating the integrals gives 



(S10) 



V^CNccl / V^Nccl + ^ Nccl ~ 7Nccl 

This is plotted for different ^Neei in inset of Fig. 2 of the main text 
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^r/YSDW 

Figure SI: Gap function of the electron pocket Ax,n as function of Matsubara frequency Lu n (in units of 7sdw) for the 
parameters used in the main text and £n oc i = (red/solid line) and £n cg i = 0.5a (blue/dashed line). 

III. BERGMANN-RAINER ANALYSIS 

We here recall how to obtain the functional derivative <5Tj! _wave /<L4]\j e ei (w) via the Bergmann-Rainer approach [f|. 
The linearized gap equation for the s-wave channel can be cast as an eigenvalue problem 



rjA = K S A 



(S12) 



with K s is a matrix in Matsubara and band space and r\ the largest eigenvalue of K . T c is the temperature at which 
77 = 0. 

Changing j4n cc ](w) — > A-^ ee \(uj) + SA^ cc \(oj) and using the Hellman-Feynman theorem gives 



Using the result: 



Sri = (aV SKs A 

<L4 Nee i(a;) V / SA Nee \(uj) 



(S13) 



da 



duj- 



ojAi (u>) 



(S14) 



we obtain: 



5A Nce i (w) 



^Nccl 



<5 m „sgn(w„)/„ (^) 



l0 2 + \LJ n -U m \^ 



(S15) 



where f n (x) = xlm [ip (— n + ix) — ip (n + 1 + ix)] and ip(x) is the digamma function. 

Using the fact that rj varies smoothly with temperature and rearranging the equation gives 



uj 8T r 



T c SA^eel (w) 



1 uj 8K S 2 
A - — A 

SA Nee \ (uj) 



7J=0 



(si6) 



where the sum over Matsubara and band indices is left implicit. Here we choose to consider the functional derivative 
with respect to A(oj)/uj because at low frequencies this is the quantity which gives the scattering rate, while if A is 
concentrated at high frequencies the ljoo gives the usual BCS logarithm. 

We see that the logarithmic derivative of the transition temperature (|S16|) is determined by the frequency dependence 
of the gap function (eigenvectors) Ai jM - explicitly in the expectation value of SK S and implicitly in the temperature 
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dependence of 77. In figure [Si] we plot the gap Ax, n of the electron pocket for the cases ^Neei = (red/solid line) 
and £ncc1 = 0.5a (blue/dashed line). The hole pocket gap displays a similar behavior. Clearly, the s-wave gap in the 
presence of Neel fluctuations decreases much faster as function of frequency than the gap of the "pure" s H state. 
Thus, the state adapts to the presence of (tt, tt) fluctuations by suppressing the gap at higher frequencies to avoid 
the main depairing effect, which comes from the high frequency components of the Neel spectrum. 



IV. DEPENDENCE OF TRANSITION TEMPERATURE ON IMPURITY SCATTERING 



The Bergmann-Rainer approach is also used to obtain the dependence of the transition temperature on impurity 
scattering: 



dr- 



df 



(S17) 



r?=0 



where t 1 is the scattering rate due to impurity scattering. The reduced wave is then obtained via the linear 
approximation T[ 



/C rj^S — V 



-1 I 9T < 



dr- 1 
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